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Abstract 

The universal behaviour of the short-time dynamics of the three state 
Potts model in two dimensions at criticality is investigated with Monte 
Carlo methods. The initial increase of the order is observed. The new 
dynamic exponent 9 as well as exponent z and (5/v are determined. The 
measurements are carried out in the very beginning of the time evolution. 
The spatial correlation length is found to be very short compared with the 
lattice size. 
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for years it has been known that there exist universahty and scahng behaviour 
for statistical systems at criticahty in equihbrium or near equihbrium, more or 
less due to the infinite spatial and time correlation length. Recently it has been 
observed that universality and scaling may also be present far from equilibrium. 
One of the examples is that the Ising model, initially in very high temperature, is 
suddenly quenched to the critical temperature and then evolves with the dynamics 
of model A. According to an argument of Janssen et al. with two-loop e-expansion 
[0, besides the well-known long-time universal behaviour, there exists another 
universal stage in the earlier time, termed critical initial slip, which sets in right 
after the microscopic time scale. The characteristic behaviour of such process is 
that, if a non-zero but small initial magnetization mo is generated in the system, 
the anomalous dimension of the operator mo gives rise to the critical increase of 
the magnetization 

M(t)~mot^ (1) 

with 9 being a new dynamic critical exponent. Detailed scaling analysis reveals [0 
that the characteristic time scale for the critical initial slip is to ~ mo^^^°, where 
Xq is the scaling dimension of mo, and related to 9 hj Xq = 9z + Pu. Interesting 
enough, it was pointed out that the exponents (3, v and z should be valued the 
same as those in the equilibrium or long-time stage of the relaxation. Previously 
9 has been measured with Monte Carlo simulation in two dimensions somehow 
indirectly from the power law decay of the autocorrelation p, and recently in 
three diemensions directly from the power law increase of the magnetization ||^. 
They are in good agreement with the result from the e-expansion. Furthermore, 
based on the scaling relation in the initial stage of the time evolution, a new 
promising way for measuring the exponents z, (3 and v from the finite size scaling 
has been proposed This indicates a possible broad application of the short- 
time dynamics. Therefore more and deeper understanding of this phenomenon 
becomes urgent. 

As far as we know, even though analytical perturbative calculations for the 
critical initial slip can be extended to the 0{N) vector model [|l], the numerical 
results are so far limited to the Ising model P, H, ^ |, |^. The purpose of 
this letter is to report systematic Monte Carlo simulations of the short-time 
behaviour of the two dimensional Potts model at criticality relaxing from high 
temperature states. A refined measurement of the exponent 9, z and jS/u from 
the power law behaviour of the physical observables is presented. It relates the 
above mentioned indirect and direct measurements of 9 each other and provides 
a consistent test of the scaling relation for the Potts model. Finite size effects 
and the week dependence of the measurement of 9 on mo are discussed and the 
spatial correlation length is computed. 

From the scaling analysis, the autocorrelation has the initial behaviour 

A{t) ~ (2) 
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Figure 1: Magnetization M(t)/M(0) vs. time in double-log scale for the initial 
value mo = 0.08. 

in case of tuq = 0. Most of the previous measurements of 6 were based on this 
relation. A disadvantage is that one has to input z to obtain 6. Since z is one 
order of magnitude bigger than 6, a small relative error of z will induce a big error 
in 6. This becomes more severe when 6 is getting smaller. A direct measurement 
from the power law increase of the order parameter can improve this situation ||^. 
As we will see later, for the Potts model we did observe the same initial increase 
of the order shown in The week dependence of the practical measurement 
of 6 on mo appears, however, more visible in the Potts model than in the Ising 
model. Therefore in this letter we consider a correction of 6 for finite mo. Since 6 
for the Potts model is smaller than that for the 2D Ising model, its measurement 
from the auto-correlation is harder. 

On the other hand, traditionally the dynamic exponent z is defined and mea- 
sured from the long-time exponential decay of the time correlation or the mag- 
netization of the systems @, ^ . Due to the critical slowing down this is somehow 
difficult. From the above discussion, however, it is easy to realize that with 6 
in hand, one can obtain z quite rigorously from the power law decay (H) of the 
autocorrelation. Since the measurement is carried out at the beginning of the 
time evolution, it is efficient. This is an alternative way to measure z from the 
short-time behaviour of the system. Compared with the method proposed in 
PI, the advantage is that the measurement can be carried out in a single lattice 
rather than by comparing two lattices. Finally the static exponent 13 /u can be 
obtained from the power law increase of the second moment 



Since 2/3/i/ is one order of magnitude smaller than z, its measurement is quite 
sensitive to the error of z. 

Figure 2: Autocorrelation A{t) vs. time in double-log scale for L = 9, 18, 36, 
72, 144 and 288. 

For the study of the short-time dynamics, usually quite a big lattice is used. 
In this letter, the finite size effects will be discussed. It turns out that for the 
measurement of 6 the finite size effect is not so big. Furthermore, the spatial 
correlation is measured and found to be very small compared with the lattice size. 
This indicates that the mechanism for the universality and scaling in short-time 
dynamics should be different from that in the equilibrium or near equilibrium. 

The Hamiltonian for the q state Potts model 



(3) 



ffi = 1, 



(4) 



<ij> 
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Figure 3: Second moment M^'^^t) / M^'^^l) vs. time in double-log scale. 



with < ij > representing nearest neighbors. In this letter we only consider the 
three state case. It is well known that for the three state Potts model the critical 
point locates at Jc = log(l + ^/3). As in case of the Ising model, initially the 
Potts model is prepared to be in a random state with a sharp magnetization mo. 
Then it is released to the evolution with the heat-bath algorithm at the critical 
temperature. We measure the time evolution of the magnetization, the second 
moment, the auto-correlation and the spatial correlation respectively 




where the average is taken over the independent intial configurations. Except 
for the magentization M{t), the above definitions are restricted here to the case 
of mo = 0. In spite of the lack of the analytic analysis, we assume that all the 
scaling properties including the increase of the order for the Ising model are valid 
also for the Potts model and test them by numerical simulation. In the same 
time the related critical exponents will be determined. 

In Fig. 1, as an example, the time evolution of the magnetization with the 
initial value mo = 0.08 for different lattice size L is displayed in a double log- 
scale to present the power law increase. It is remarkable that the power law 
increase starts also from the very beginning of the time evolution t = 1 as it is 
in the three dimensional Ising model. 9 can be estimated from the slope of the 
curves. It is clearly seen that 9 converges to a definite value when the lattice 
size L > 36. In other words, for the measurement of 9 the finite size effect is 
already quite small for a lattice size L — 36. In comparison to this, to observe 
the power law decay for the auto-correlation one needs much bigger lattices, as 
will be seen later. In Tab. 1, 9 for L = 72 measured from t = 1 to t = 15 
for different initial magnetizations has been summarized. The total number of 
samples for the independent initial configurations is 80, 000 for bigger mo and 
480, 000 for smaller mo. The errors are estimated by dividing the data into four 
or six groups, respectively. Different from the Ising model, the measured 9 shows 
slightly more dependence on mo. Therefore, according to its definition, a linear 
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mo 


0.10 


0.08 


0.06 


0.04 


0.00 


e 


0.1076(08) 


0.1036(12) 


0.0980(06) 


0.0925(14) 


0.0815(27) 



Table 1: The measured 9 for L = 12. 

extrapolation of 6 to the fixed point itlq = is carried out. This leads to the 
value 

e = 0.0815(27). 



L 


72 


144 


288 


oo 


-d/z + e 

{d-2(3/u)lz 


-0.8510(08) 
0.7921(11) 


-0.8387(09) 
0.7881(14) 


-0.8335(09) 
0.7875(28) 


-0.8283(20) 
0.7878(16) 



Table 2: Results for L = 72, 144 and 288. 



Now we set mo = and measure the auto-correlation. In Fig. 2, the depen- 
dence of the auto-correlation on L is presented. Obviously at L = 36 no power 
law behaviour is observed. The convergence to a power law behaviour only starts 
around L = 144. It is clear that the regime presenting power law grows when 
the lattice size increases. It is interesting that, somehow different from M{t), the 
first time steps apparently deviate slightly from the power law. In Tab. 2 the 
corresponding values for —d/z + 6 measured from t = 5 to t = 50 are given. We 
stop the measurement at t = 50 due to the obvious bigger finite size effect and 
statistical errors. The total samples for L = 144 is 40, 000 and for L = 288 is 
16, 000. if we only intend to obtain z by taking 6 as input, we should already 
be satisfied with the lattice size since the results from L = 144 and L = 288 are 
already so close. However, in order to get the better /3/z/ later by inputting the 
z measued here, we perform the linear extrapolation for —d/z + 6 over 1/L to 
L = oo and obtain 

z = 2.1983(81). 

Compared to the values of z distributed between z = 2.2 and z = 2.7 from 
different numerical measurments [|ll], 0, [l^, |14| , our result supports the relative 



small z iTTI, n!2 



In Fig. 3, the power law increase of the second moment M^'^'>(t) is shown. The 
measurement of {d — 2j3/h')/z from t = 5 to t = 50 are also listed in Tab. 2. From 
the data we can see that for L = 144 and L = 288, the finite size effects are 
already less prominent than the statistical errors. Therefore the value 0.7878(22) 
of [d — 2(3 / u) / z a.t L = oo is only a simple average of them. Here we get 



2l3/u = 0.2682(73), 
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which is in good agreement with the exact value 4/15 ~ 0.2667 |T0[. Such 
coincidence provides a strong support for scahng in the short-time dynamics. 

Finally we have also measured the correlation length ^(t) from the spatial 
correlation function C{x,t). for example, for L = 288, ^(t = 96) ~ 6.0 and it is 
much smaller than the lattice size L. 

In conclusion, by observing the power law behaviour of the magnetization 
M(t), the second moment M^'^\t) and the auto-correlation A{t), we confidently 
confirm the scaling properties for the Potts model in the short-time dynamics and 
obtain the related critical exponents 6, z and /3/z/. This is the first measurment of 
Q for the Potts model. Our way to determine z is efficient. Such an investigation 
for models in other universality classes should be carried out. 
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